Cementation of kerogen-rich marls by alkaline fluids released during weathering of thermally metamorphosed marly sediments. Part II: Organic matter evolution,magnetic susceptibility and metals (Ti,Cr, Fe) at the Khushaym Matruk natural analogue (Central Jordan)

Spontaneous combustion, less than 1 Ma ago, affected a 60-m thick sediment pile of biomicrite at the Khushaym Matruck site (Jordan). The present study shows that three retrograde alteration stages occurred: weathering, thermal stress and oxidative alkaline perturbation. μ-FT-i.r. spectra of isolated kerogens and oxygen index of whole rocks indicate that oxidation of organic matter occurred down to ∼10 m beneath the metamorphosed zone at Khushaym Matruck. The occurrence of the oxidative weathering bacterially mediated, as suggested by the mass chromatograms of saturated hydrocarbons, can explain high Rock-Eval T_max values and low petroliferous potential measured along the sedimentary pile. On the other hand, the thermal extent of combustion events was limited to the first 2 m from the contact. The mean reflectance of 0.20–0.24% and porosity of ca. 50% of the grey clayey biomicrites indicate that organic matter was very immature and sediments were unconsolidated at the time of the combustion event. Using mineralogy, microscopic analyses of vegetable debris and magnetic susceptibility, a suite of characteristic points corresponding to the thermal imprint can be assessed: (i) x = 0m, T ∼ 1000 °C, (ii) x = 1 m, T ∼ 350 °C, (iii) x = 2 m, T ∼ 150 °C and (iv) x > ∼ 8 m, T ∼ 30 °C. Paleocirculation of meteoric groundwater in the ‘cement-marbles’ generated high-pH fluids that have circulated via fractures and through the matrix porosity of the underlying biomicrites but have also induced alkaline hydrolysis and oxidative attack of the organic matter. The polysaccharide/lignin ratio derived from μ-FT-i.r. analyses shows that the delignification of vegetable debris and degradation of polysaccharides progressively decline in the indurated zone, which indicates a decrease in the pH of migrating solutions. The latter also severely oxidized organic matter at 2.10 and 3.05 m as revealed by the oxygen index and induced the generation of bitumen. The spatial correlation between the oxidation levels of organic matter and the metal contents (Fe, Ti and Cr) suggests that redox reactions were responsible for the immobilization of metals in the indurated biomicrites. The intensity of these reactions is attributed to changes in the fluid flow regime within the sedimentary column.     

An experimental study of Ostwald ripening of olivine and plagioclase in silicate melts: implications for the growth and size of crystals in magmas

We carried out an experimental study to characterize the kinetics of Ostwald ripening in the forsterite-basalt system and in the plagioclase (An₆₅)-andesite system. Eight experiments were done in each system to monitor the evolution of mean grain size and crystal size distribution (CSD) with time t; the experiments were performed in a 1-atmosphere quench furnace, at 1,250°C for plagioclase and 1,300°C for olivine. Very contrasted coarsening kinetics were observed in the two series. In the plagioclase series, the mean grain size increased as log(t), from ≈3 μm to only 8.7 μm in 336 h. The kinetic law in log(t) means that Ostwald ripening was rate-limited by surface nucleation at plagioclase-liquid interfaces. In the olivine series, the mean grain size increased as t ^1/3, from ≈3 μm to 23.2 μm in 496 h. A kinetic law in t ^1/3 is expected when Ostwald ripening is rate-limited either by diffusion in the liquid or by grain growth/dissolution controlled by a screw dislocation mechanism. The shape of olivine CSDs, in particular their positive skewness, indicates that grain coarsening in the olivine experiments was controlled by a screw dislocation mechanism, not by diffusion. As the degrees of undercooling ΔT (or supersaturation) involved in Ostwald ripening are essentially <1°C, the mechanisms of crystal growth identified in our experiments are expected to be those prevailing during the slow crystallisation of large magma chambers. We extrapolated our experimental data to geological time scales to estimate the effect of Ostwald ripening on the size of crystals in magmas. In the case of plagioclase, Ostwald ripening is only efficient for mean grain sizes of a few microns to 20 μm, even for a time scale of 10⁵ years. It can, however, result in a significant decrease of the number of small crystals per unit volume, and contribute to the development of convex upwards CSDs. For olivine, the mean grain size increases from 2–3 μm to ≈70 μm in 1 year and 700 μm in 10³ years; a mean grain size of 3 mm is reached in 10⁵ years. Accordingly, the rate of grain size-dependent processes, such as compaction of olivine-rich cumulates or melt extraction from partially molten peridotites, may significantly be enhanced by textural coarsening.     

Projection of subtropical gyre circulation and associated sea level changes in the Pacific based on CMIP3 climate models

For all of the IPCC Special Report on Emission Scenarios (SRESs), sea level is projected to rise globally. However, sea level changes are not expected to be geographically uniform, with many regions departing significantly from the global average. Some of regional distributions of sea level changes can be explained by projected changes of ocean density and dynamics. In this study, with 11 available Coupled Model Intercomparison Project Phase 3 climate models under the SRES A1B, we identify an asymmetric feature (not recognised in previous studies) of projected subtropical gyre circulation changes and associated sea level changes between the North and South Pacific, through analysing projected changes of ocean dynamic height (with reference to 2,000 db), depth integrated steric height, Sverdrup stream function, surface wind stress and its curl. Poleward expansion of the subtropical gyres is projected in the upper ocean for both North and South Pacific. Contrastingly, the subtropical gyre circulation is projected to spin down by about 20 % in the subsurface North Pacific from the main thermocline around 400 m to at least 2,000 m, while the South Pacific subtropical gyre is projected to strengthen by about 25 % and expand poleward in the subsurface to at least 2,000 m. This asymmetrical distribution of the projected subtropical gyre circulation changes is directly related to differences in projected changes of temperature and salinity between the North and South Pacific, forced by surface heat and freshwater fluxes, and surface wind stress changes.     

Report on the 4th International Meeting of the IUGS Lower Cretaceous Ammonite Working Group, the “Kilian Group” (Dijon, France, 30th August 2010)

The 4th Kilian Group meeting (Dijon, France, 30th August 2010) focused on the Aptian and Albian Stages. For the Aptian, a two-fold division of the stage was adopted for the Mediterranean area with a boundary between the Dufrenoyia furcata and Epicheloniceras martini Zones. The main changes to the zonal scheme concern the Lower Aptian with: the introduction of a Deshayesites luppovi Subzone in the upper part of the Deshayesites oglanlensis Zone; the replacement of Deshayesites weissi by Deshayesites forbesi as new index-species of the second interval zone; the introduction of a Roloboceras hambrovi Subzone in the upper part of the D. forbesi Zone; and the subdivision of the D. furcata Zone into the D. furcata and Dufrenoyia dufrenoyi Subzones. For the Albian, the upper part of the Douvilleiceras mammillatum Zone (Lower Albian) is now characterized by a Lyelliceras pseudolyelli Subzone. The main amendments concern the Upper Albian. The base of this substage is defined by the base of the Dipoloceras cristatum Zone. Above it, the Upper Albian zonal scheme comprises in stratigraphic order the Mortoniceras pricei, Mortoniceras inflatum, Mortoniceras fallax, Mortoniceras rostratum, Mortoniceras perinflatum and Arrhaphoceras briacensis Zones.     

Approximate Derivative Computations for the Gradient-Based Optimization Methods in the Local Gradual Deformation for History Matching

Reservoir characterization needs the integration of various data through history matching, especially dynamic information such as production or four-dimensional seismic data. To update geostatistical realizations, the local gradual deformation method can be used. However, history matching is a complex inverse problem, and the computational effort in terms of the number of reservoir simulations required in the optimization procedure increases with the number of matching parameters. History matching large fields with a large number of parameters has been an ongoing challenge in reservoir simulation. This paper presents a new technique to improve history matching with the local gradual deformation method using the gradient-based optimizations. The new approach is based on the approximate derivative calculations using the partial separability of the objective function. The objective function is first split into local components, and only the most influential parameters in each component are used for the derivative computation. A perturbation design is then proposed to simultaneously compute all the derivatives with only a few simulations. This new technique makes history matching using the local gradual deformation method with large numbers of parameters tractable.     

Initialisation and predictability of the AMOC over the last 50 years in a climate model

The mechanisms involved in Atlantic meridional overturning circulation (AMOC) decadal variability and predictability over the last 50 years are analysed in the IPSL–CM5A–LR model using historical and initialised simulations. The initialisation procedure only uses nudging towards sea surface temperature anomalies with a physically based restoring coefficient. When compared to two independent AMOC reconstructions, both the historical and nudged ensemble simulations exhibit skill at reproducing AMOC variations from 1977 onwards, and in particular two maxima occurring respectively around 1978 and 1997. We argue that one source of skill is related to the large Mount Agung volcanic eruption starting in 1963, which reset an internal 20-year variability cycle in the North Atlantic in the model. This cycle involves the East Greenland Current intensity, and advection of active tracers along the subpolar gyre, which leads to an AMOC maximum around 15 years after the Mount Agung eruption. The 1997 maximum occurs approximately 20 years after the former one. The nudged simulations better reproduce this second maximum than the historical simulations. This is due to the initialisation of a cooling of the convection sites in the 1980s under the effect of a persistent North Atlantic oscillation (NAO) positive phase, a feature not captured in the historical simulations. Hence we argue that the 20-year cycle excited by the 1963 Mount Agung eruption together with the NAO forcing both contributed to the 1990s AMOC maximum. These results support the existence of a 20-year cycle in the North Atlantic in the observations. Hindcasts following the CMIP5 protocol are launched from a nudged simulation every 5 years for the 1960–2005 period. They exhibit significant correlation skill score as compared to an independent reconstruction of the AMOC from 4-year lead-time average. This encouraging result is accompanied by increased correlation skills in reproducing the observed 2-m air temperature in the bordering regions of the North Atlantic as compared to non-initialized simulations. To a lesser extent, predicted precipitation tends to correlate with the nudged simulation in the tropical Atlantic. We argue that this skill is due to the initialisation and predictability of the AMOC in the present prediction system. The mechanisms evidenced here support the idea of volcanic eruptions as a pacemaker for internal variability of the AMOC. Together with the existence of a 20-year cycle in the North Atlantic they propose a novel and complementary explanation for the AMOC variations over the last 50 years.     

Decadal fingerprints of freshwater discharge around Greenland in a multi-model ensemble

The recent increase in the rate of the Greenland ice sheet melting has raised with urgency the question of the impact of such a melting on the climate. As former model projections, based on a coarse representation of the melting, show very different sensitivity to this melting, it seems necessary to consider a multi-model ensemble to tackle this question. Here we use five coupled climate models and one ocean-only model to evaluate the impact of 0.1 Sv (1 Sv = 10⁶ m³/s) of freshwater equally distributed around the coast of Greenland during the historical era 1965–2004. The ocean-only model helps to discriminate between oceanic and coupled responses. In this idealized framework, we find similar fingerprints in the fourth decade of hosing among the models, with a general weakening of the Atlantic Meridional Overturning Circulation (AMOC). Initially, the additional freshwater spreads along the main currents of the subpolar gyre. Part of the anomaly crosses the Atlantic eastward and enters into the Canary Current constituting a freshwater leakage tapping the subpolar gyre system. As a consequence, we show that the AMOC weakening is smaller if the leakage is larger. We argue that the magnitude of the freshwater leakage is related to the asymmetry between the subpolar-subtropical gyres in the control simulations, which may ultimately be a primary cause for the diversity of AMOC responses to the hosing in the multi-model ensemble. Another important fingerprint concerns a warming in the Nordic Seas in response to the re-emergence of Atlantic subsurface waters capped by the freshwater in the subpolar gyre. This subsurface heat anomaly reaches the Arctic where it emerges and induces a positive upper ocean salinity anomaly by introducing more Atlantic waters. We found similar climatic impacts in all the coupled ocean–atmosphere models with an atmospheric cooling of the North Atlantic except in the region around the Nordic Seas and a slight warming south of the equator in the Atlantic. This meridional gradient of temperature is associated with a southward shift of the tropical rains. The free surface models also show similar sea-level fingerprints notably with a comma-shape of high sea-level rise following the Canary Current.